Ensemble averages and nonextensivity at the edge of chaos of one-dimensional maps 
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Ensemble averages of the sensitivity to initial conditions ^(t) and the entropy production per 
unit time of a new family of one-dimensional dissipative maps, Xt+i = 1 — ae"^'''^'' {z > 0), 
and of the known logistic-like maps, xt+i = 1 — a\xt\^{z > 1), are numerically studied, both for 
strong (Lyapunov exponent Ai > 0) and weak (chaos threshold, i.e., Ai — 0) chaotic cases. In 
all cases we verify that (i) both (In, ^) [In^ a:: = (a::^"' — 1)/(1 — g); lui a; = Ina;] and (Sg) [Sq = 
(1 — pDI (q — 1); S\ = — pi Inpi] linearly increase with time for (and only for) a special value 
of q, qsen, and (ii) the slope of (In^ ^) and that of (Sq) coincide, thus interestingly extending the well 
known Pesin theorem. For strong chaos, q^en = Ij whereas at the edge of chaos, qgen{z) < 1- 
PACS number: 05.20.-y, 05.45.-a, 05.45.Ac, 05.45. Df 



Low-dimensional non-linear maps play an important 
role in the development of the theory of chaos in physics 
and mathematics. They exhibit various routes to chaos 
and their related metric universality classes [1]. In par- 
ticular, one-dimensional maps are paradigmatic models 
to study the emergence of complexity in dynamical sys- 
tems. In this field, much work has been done to find 
properties of chaos that could enable a classification of 
deterministic systems. In particular, dynamical indica- 
tors like sensitivity to initial conditions, Lyapunov expo- 
nents, Kolmogorov-Sinai (KS) entropy [2,3], topological 
entropy and others have been developed. Nevertheless, 
it is known today that there are ubiquitous natural and 
artificial dynamical systems for which these standard in- 
dicators give but a poor description of the complexity 
of their time evolution. This is typically the case in the 
frontier between standard chaos and regular orbits. Such 
complex deterministic behaviour has been termed weak 
chaos, to distinguish it from strong or full chaos. The 
latter correspond to positive metric entropy or, equiva- 
lently, positive Lyapunov exponent. The former occurs 
for zero values of these indicators. This phenomenon has 
led to the generalization of the involved concepts [4]. 

A new and fruitful approach that characterizes these 
weakly chaotic dynamical systems is related to the nonex- 
tensive generalization [5,6] of the Boltzmann-Gibbs (BG) 
statistical mechanics. This generalization has since more 
than one decade raised much interest in physical situa- 
tions that do not satisfy the customary BG thermal equi- 
librium conditions, e.g., quick mixing in phase space and 
ergodicity. In such cases, anomalous dynamical proper- 
ties are the rule: see [7-9], among many others. 

In the present paper, we study numerically the sen- 
sitivity to initial conditions and the entropy production 
per unit time (a concept completely analogous to the 
Kolmogorov-Sinai entropy rate) of the standard z-logistic 
maps as well as of a new family of one-dimensional dis- 
sipative maps, from now on referred as to the exponen- 



tial maps, here introduced in order to have one more 
paradigmatic universality class on which basic proper- 
ties can be tested. In all the calculations we present 
here, we make ensemble averages by sampling randomly 
the entire phase space. This procedure will be shown to 
substantially modify the behavior at the edge of chaos, 
whereas nothing very drastic is changed for strong chaos. 
These results lead us to the numerical confirmation, for 
both families of maps, of an appropriate generalization 
of the well known Pesin identity (i.e., the coincidence of 
the positive Lyapunov exponent and KS entropy rate). 
This generalization was conjectured in 1997 [4], recently 
proved for special trajectories (more precisely as upper 
bounds of all trajectories; see [10] for the logistic maps), 
and is here verified, for the first time, for ensemble aver- 
ages. 

Let us first recall the z-logistic family of maps 

xt+i = I - a\xt\'' (z > l;ae [0,2];|a;t| < 1). (1) 

The critical value adz) {chaos threshold or edge of chaos) 
monotonically increases from 1 to 2 when z increases from 
1 to infinity. The z = 2 map is, as well known, isomor- 
phic to yt+i oc j/f (1 — yt). The z-logistic maps constitute 
important universality classes of unimodal maps. The 
larger z is, the flatter is the a; = maximum. To char- 
acterize a degree of flatness that z ^ oo cannot attain, 
we introduce the following new family of maps, inspired 
in Cauchy's exponential function (infinitely differentiable 
at a; = and nevertheless nonanalytic): 

xt+i = 1 - ae~i/l"'l^ (z > 0; a e [0, a*{z)]; \xt\ < 1) , (2) 

where a*{z) depends slowly from z (e.g., a* (0.5) ~ 5.43). 
We only address here z > Zc — 0.4, Zc being the value 
above which the attractors are topologically isomorphic 
(see Fig. 1) to those of the logistic map. In fact, one 
observes a whole cascade of period doublings 2*^^^ 2^ 
at parameter values Uk that accumulate at a^z) (e.g., 
ac(0.5) = 3.32169594...) above which chaos exists. 
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FIG. 1. a-dependence of the dynamical attractor of the 
z = 0.5 exponential map. 



Table- Logistic (top) and exponential (bottom) maps. 



z 




av 
^sen 






1.10 


1.12498... 


0.43 ±0.01 


0.12 ±0.01 


0.12 ±0.02 


1.25 


1.20951... 


0.42 ± 0.01 


0.18 ±0.01 


0.18 ±0.02 


1.50 


1.29550... 


0.40 ±0.01 


0.22 ±0.01 


0.22 ±0.02 


1.75 


1.35506... 


0.37 ±0.01 


0.26 ±0.01 


0.26 ±0.02 


2.00 


1.40115... 


0.36 ±0.01 


0.27 ±0.01 


0.27 ±0.02 


2.50 


1.47055... 


0.34 ±0.01 


0.28 ±0.01 


0.28 ±0.02 


3.00 


1.52187... 


0.32 ±0.01 


0.29 ±0.02 


0.29 ±0.03 


5.00 


1.64553... 


0.28 ±0.01 


0.30 ±0.02 


0.30 ±0.03 


0.40 


3.05996... 


0.40 ±0.01 


0.25 ±0.01 


0.25 ±0.02 


0.50 


3.32169... 


0.35 ±0.01 


0.27 ±0.01 


0.27 ±0.02 


0.75 


3.68229... 


0.30 ± 0.01 


0.29 ± 0.02 


0.29 ±0.02 


1.00 


3.9()7()3... 


0.23 ± 0.02 


0.32 ± 0.02 


0.32 ± 0.02 


1.25 


4.07088... 


0.20 ±0.02 


0.39 ±0.02 


0.39 ±0.03 


1.50 


4.19820... 


0.15 ±0.02 


0.43 ± 0.02 


0.43 ±0.03 



The ^;-logistic maps, and many others, have al- 
ready been deeply studied in the regions of strong 
chaos. The sensitivity to the initial conditions £,{t) = 
l™Arr(o)^o Ax\o) ' ^^crc Ax(0) is the discrepancy of ini- 
tial conditions at time t = and Ax(t) its time de- 
pendence, satisfies the differential equation d^/dt = Ai^, 
where Ai is the Lyapunov exponent, hence ^ = e^^*. Con- 
sequently, if Ai < (Ai > 0), the system is said to be 
strongly insensitive (sensitive) to the initial conditions 
(and intermediate rounding). If, however. Ai = (as it 
happens at the edge of chaos), then ^(<) is expected to 
satisfy a more general equation, namely d^/dt = Ag^' (Xq 
being a (/-generalized Lyapunov coefficient), hence [4] 

^=[l + {l-q)Xgt]^ =e^^K (3) 

This function, q- exponential, recovers the exponential 
function (a power law) for g = 1 (g ^ 1.) If g < 1 



and Ag > (g > 1 and Ag < 0), the system is said to be 
weakly sensitive [insensitive) to the initial conditions. 

Consistently with these results, whenever strong chaos 
is present, the appropriate entropy is known to be the BG 
one, Sbg = ^ Sl^ift InPi- But, at the edge of chaos, 
Sbg becomes poorly informative, and we must instead 
use ( [4,5,11-13] and references therein) 

Sg= iS^ = SBG)- (4) 

This nonextensive entropy enables in fact the generaliza- 
tion [5] of BG statistical mechanics itself. For each uni- 
versality class of edge of chaos, one special value of the en- 
tropic index q is to be used, which we note qsen, where sen 
stands for sensitivity. Up to now, four different methods 
have been developed which provide the special value qsen- 
The first method is based on the sensitivity to initial con- 
ditions [4,11], the second one is related to the entropy 
production per unit time [14,15], the third one is based on 
the geometrical description of the multifractal attractor 
[12], and the last one concerns relaxation properties [16]. 
We address here the first two methods. For the 2:-logistic 
map, qsen{z) is known for typical values of z [11,12]: Sec 
filled circles in Fig. 3 (e.g., qsen{z) = 0.2445...; also [10] 
Ao.2445 = In a/ In 2 = 1.3236..., where a is the Feigen- 
baum constant). But. for the present exponential family, 
this is still unpublished and out of the scope of the present 
work. Our main goal here is to study what happens when 
we do averages (noted (...)) over the entire phase space 
(i.e., — 1 < a; < 1), instead of using the special region 
around a; = (location of the maximum of the map) . 

For the sensitivity we proceed as follows. We consider, 
at an initial location of X, two points very close (say at a 
distance Aa;(0) = 10~^^) and then numerically calculate 
^(t) from its definition. We do this operation many times, 
starting from values of x randomly chosen in the allowed 
interval. Finally, we average all the values of lug ^, where 
the inverse function of the g-exponential function is de- 
fined as follows: 

- 1 

lUg X = — (Ini a; = In a;) . (5) 

Various values of q are tested until we obtain a linear 
time dependence of (lng^)(t). That is the special value 
of q we are looking for, and we denote it qg^,^, where 
av stands for average; we denote A^a^ the linear coeffi- 
cient. See Fig. 2. For example, for the edge of chaos 
of the z = 2 logistic map. we obtain g^J,^ ~ 0.36, in re- 
markable contrast with qsen — 0.2445. We shall come 
back onto this point later. No such discrepancy is ob- 
tained for strong chaos, where we verify in all cases that 
<?sera — Isen = 1- The Same method is applied to the ex- 
ponential family of maps: See Fig. 2 (e.g., at the edge of 
chaos ac(0.5) = 3.32169594, we obtain - 0.35). 
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FIG. 2. Time dependence of (In, ^) and (Sq): z = 2 logistic map for strong [(a) a = 2] and weak [(c) a = 1.401155189] chaos, 
and z = 0.5 exponential map for strong [(b) a = 4] and weak [(d) a — 3.32169594] chaos. Sensitivity function {lnq^){t): 
averages over lO'^ (10'^) runs for (a) and (6) ((c) and (d)); we use Ax(0) = 10 as the initial discrepancy unless otherwise 
indicated; in the insets, we show the Zmear tendency of the sensitivity function for q^^^ with various values of Aa;(0); at the edge 
of chaos ((c) and (d)) we exhibit the q — 1 curve nonlinearity. Entropy {Sq){t): (a,b) 3000 runs with A'^ = WW with W — 10® 
and W = (empty and filled symbols respectively); 50000 runs with N = lOW with W = 10® (for (c)) and W = 5.10* and 

W = 10® (for (d)). (c) inset: determination of g"e„ (see text), (d) inset: we exhibit the g = 1 curve nonlinearity. 
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FIG. 3. 2;-dependence of g^e^ (empty circles and squares: 
present work) and qsen (filled circles: from [11,12]; filled 
squares: from [18]). Dotted lines are guides to the eye. 

For the entropy production wc proceed as introduced 
in [17] for conservative systems and then used in [14] for 
the logistic map. We divide the phase space in W equal 
intervals, and put initially (randomly or not) points 
in one of them. We then accompany the spread of points 
within phase space, and calculate Sq{t). This enables the 
calculation of the entropy production per unit time 



Kq = lim lim lim 



(6) 



We then repeat many times starting from randomly cho- 
sen intervals and finally average the entropy Sq itself. We 
test various g''s until we obtain a linear time dependence 
of {Sq){t); we note K^av the linear coefficient. It turns 
out that in all cases wc obtained the same value q^^^^ pro- 
vided by the method of the sensitivity (see Fig. 3). And 
in all cases we obtained K^i'v = A^a„ , which remarkably 
extends the standard Pesin theorem [19]. 

This kind of averaging, on one hand mimics Gibbs' ap- 
proach of thermostatistical ensembles, on the other hand 
is a natural way to minimize the considerable fluctua- 
tions which appear at the edge of chaos a = ac{z). At 
this threshold, and in contrast with the case where strong 
chaos exists, considerably large fluctuations appear in say 
Sq{t). In order to get rid of these fluctuations, we average 
{Sq){t) over a large enough number of initial windows in- 
stead of only using the "best windows" as done in [14,15]. 
Here, we typically used W/2 windows chosen randomly. 
The fluctuations that so remain are considerably smaller. 

In the chaotic regime for the z = 2 logistic map we 
obtained, for a = 2, the well known result Ai = Ki = 
ln2 ~ 0.69. In the chaotic regime (a — 4) for the z = 0.5 
exponential map, we obtained Ai ~ iiTi ~ 0.38. 

In order to make the entropy result more precise (see 
the inset (c) of Fig. 2 (bottom) for the standard logistic 
map as illustration) we fitted the curves {Sq){t) in the 



time interval [ti, with the polynomial A + Bt + Ct^. 
We define the nonlinearity measure R = C{ti + t2)/B 
{R = for a perfect straight line). For the logistic map 
(z = 2), wc choose (ti,t2) = (1,30) for all g's, so that 
the factor {ti + t2) is just a normalizing constant. This 
inset shows that R vanishes for q = 0.36 ± 0.01. These 
results are not sensitive to changes in t2 < 30. We used 
the same procedure for the exponential maps. 

Let us summarize our main results: (i) We intro- 
duced a new universality class for one-dimensional uni- 
modal dissipative maps, Eq. (2), which corresponds to a 
proper characterization of extremely flat maps {z ^ oo 
of Eq. (1)); (ii) For this family of maps as well as 
for the 2;-logistic ones, we performed Gibbs-ensemble- 
like averages on both the sensitivity and of the entropy 
production, and determined (exclusively from dynam- 
ics) the cntropic index qg^niz)i one and the same from 
both procedures (as it is known to occur for qsen{z))', (hi) 
For strong chaos, we verify q^eni^) = qsen{z) = 1 and 
KZSz) = A-„(^) = Kq_{z) = A,„„(z); (iv) At the 
edge of chaos, we verify that (z) decreases for increas- 
ing flatness (i.e., increasing z, in contrast with qseniz)), 
and KZAz) = XZ„iz) < Kq_{z) = Ag^„(z); the fact 
that generically qseni^) <lsen{z) is a direct consequence 
of the fact that different initial conditions yield strongly 
fluctuating trajectories which follow along the path of 
virtually any other trajectory but at shifted times [10] 
(no such phenomenon exists for strong chaos). All these 
results are expected to contribute to the correct inter- 
pretation of experimental results in dynamical complex 
systems, in particular the ubiquitous dissipative ones. 
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